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I J.  ABSTRACT 

A  12-month  program  with  the  objective  to  Develop  a  mathematical  model  of 
fluid  flow  in  Coanda  eductors  is  about  50%  complete.  The  model  is  to  serve 
as  a  tool  in  design  of  Coanda  eductors  for  pneumatic  conveying  of  excavated 
rock  materials.  Based  on  work  completed,  it  has  been  concluded  that  Coanda 
eductors  possess  potential  advantages  compared  to  alternate  ejectors  or 
other  momentum  transfer  devices.  The  advantages  are  primarily  related  to 
the  simplicity  of  the  mechanical  design./'  A  review  of  literature  has  shown 
that  very  Tittle  information  is  available  onicerning  the  fluid  mechanics  of 
the  Coanda  eductor,  the  needed  information  will  have  to  be  developed  as  part 
of  the  present  program.  Of  the  alternate  mathematical  approaches  available, 
one  based  on  application  of  an  existing  fluid  mechanics  code  was  considered 
the  most  desirable,  and  continuing  work  will  be  keyed  to  this  mathematical 
approach.  Ongoing  work  includes  development  of  the  fluid  mechanics  code, 
and  experimental  measurements  of  flow  fields  in  model  eductors  which  will 
provide  input  data  for  the  model  and  serve  as  a  basis  for  evaluating 
model  predictions. 
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DEVELOPMENT  OF  A  CONCEPT  FOR  A  HIGH  CAPACITY 
PNEUMATIC  CONVEYING  SYSTEM  EMPLOYING  A  FLUID  ATTACHMENT  DEVICE 
FOR  USE  IN  UNDERGROUND  EXCAVATION 
to 

U.S.  BUREAU  OF  MIfJES 
CONTRACT  NUMBER  H0220027 
INTRODUCTION 

This  report  describes  the  work  completed  during  the  first  six  months 
of  a  program  being  carried  out  by  Battelle-Northwest  for  the  U.S.  Bureau 
of  Hines.  The  objective  of  this  research  program  is  to  develop  a  mathe¬ 
matical  model  to  describe  the  fluid  mechanics  of  a  Coanda  eductor.  A 
satisfactory  model  will  serve  as  a  tool  In  design  of  Coanda  eductors  for 
pneumatic  conveying  of  excavated  rock  materials.  The  report  includes  a  sue 
aa-y  of  work  accomplished,  conclusions  and  recommendations,  followed  by 
detailed  discussions  presented  in  sections  entitled  Literature  Review, 
Mathematical  Model  Development  and  Experimental  Program. 

The  possible  use  of  a  Coanda  eductor  for  pneumatic  transport  was  sug¬ 
gested  by  preliminary  experiments  carried  out  at  Battelle-Northwest  A 
small  model  was  found  to  be  highly  effective  In  entraining  granular  mate¬ 
rials  and  discussions  with  personnel  from  the  U.S.  Bureau  of  Mines  led 
to  the  current  Interest  In  assessing  the  potential  use  of  such  an  eductor 
In  transporting  rock  in  underground  excavations. 

The  term  "Coanda"  derives  from  the  Romanian-born  engineer  named 
Henri  Coanda.  Henri  Coanda  nearly  suffered  an  early  demise  as  a  result 
of  unexpected  flow  along  exhaust  deflector  plates  he  had  Installed  on  a 
wooden  airplane  powered  by  a  type  of  jet  engine.  Instead  of  deflecting 
the  exhaust  gas  away  from  the  wooden  fuselage.  Installed  deflector  plates 
were  actually  sucking  the  flames  toward  It.  Coanda  survived  this  mishap, 
and  spent  much  of  the  remainder  of  his  life  studying  the  curious  attach¬ 
ment  phenomena  which  accompany  flow  along  solid  surfaces.  Most  of  the 
scientific  effort  devoted  to  attached  or  “Coanda"  flow  has  been  directed 
to  the  use  of  the  phenomenon  In  vehicle  propulsion. 
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Use  of  the  fluid  attachment  principle  in  an  eductor  is  attractive 
because  of  very  simple  mechanical  designs  which  are  possible,  and  more 
importantly,  because  of  the  possibility  of  designing  eductors  with  large 
throat  openings  free  from  obstructions.  These  advantages  are  of  poten¬ 
tially  great  Importance  In  pneumatic  conveying  of  solids.  This  problem 
had  not  been  studied  previously  to  an  appreciable  extent.  The  approach  taken 
In  the  present  study  is  to  develop  a  mathematical  model  which  adequately 
describes  fluid  flow  in  Coanda  eductors.  The  merits  for  using  a  Coanda 
eductor  in  a  specific  case  can  then  be  judged  from  eductors  designed  to 
meet  specific  requirements  of  flow  rate  and  pressure  drop. 

SUMMARY 


A  12-month  program  with  the  objective  to  develop  a  mathematical  model 
of  fluid  flow  In  Coanda  eductors  is  about  50%  complete.  Tasks  which  have 
been  completed  include  writing  a  detailed  program  plan,  reviewing  available 
technical  literature,  and  making  a  comparison  of  performance  of  Coanda  eductors 
with  other  momentum  transfer  methods.  Ongoing  work  includes  development 
of  a  numerical  computer  program  for  solving  the  differential  equation  of 
fluid  flow  for  axial  symmetric  flow,  and  experimental  measurements  of  flow 
fields  In  model  eductors. 

Many  studies  of  wall  jet  flow  have  been  performed  for  two-dimensional 
flow  over  flat  and  curved  surfaces.  The  experimental  and  theoretical  data 
that  are  available,  however,  do  not  apply  directly  to  flow  In  cylindrical, 
axl symmetric  geometries  of  Interest  for  the  Coanda  eductor. 

Comparison  of  the  Coanda  eductor  with  other  types  of  momentum  trans¬ 
fer  devices  indicated  that  the  chief  advantages  offered  by  this  eductor 
lie  In  Its  extreme  mechanical  simplicity,  and  In  Its  large  open  throat. 

Energy  requirements  for  the  Coanda  eductor  are  expected  to  be  higher  than 
for  mechanical  systems  such  as  conveyers  or  airlock  feeders.  Energy 
requirements  are  expected  to  be  equivalent  to  those  for  presently  avail¬ 
able  ejectors. 
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Three  approaches  to  building  a  mathematical  model  for  the  flow  field 
were  considered.  These  included  model  theory  based  on  dimensional  analysis, 
similarity  theory,  and  numerical  solution  of  the  complete  Navier-Stoakes 
equations.  A  decision  was  made  to  pursue  the  latter  approach  because  it 
can  provide  the  most  detailed  description  of  the  flow  field  with  the 
least  amount  of  experimental  data.  An  existing  computer  code  which  solves 
the  Navier-Stokes  equations  for  flows  having  axial  symmetry  was  modified 
to  calculate  the  flow  field  for  a  Coanda  eductor.  To  the  present,  steady- 
state,  well  converged  solutions  have  been  obtained  using  a  computational 
grid  with  1200  nodes.  Effort  in  the  immediate  future  will  concentrate 
on  solution  of  the  pressure  equation  and  on  means  for  realistically 
handling  eddy  viscosities. 

An  experimental  apparatus  has  been  designed  to  provide  data  on  pres¬ 
sure  drop,  entrained  flow,  and  velocity  profiles  for  several  geometries. 

The  experimental  data  will  be  used  to  test  the  validity  of  the  mathemati¬ 
cal  model  and  to  provide  Input  regarding  numerical  values  of  eddy 
viscosity. 


rnNCLUSIONS  AND  RECOMMENDATIONS 

The  work  completed  to  date  supports  the  following  conclusions: 

.  The  Coanda  eductors  possess  potential  advantages  compared  to  alternate 
ejectors  or  other  momentum  transfer  devices.  The  advantages  are 
primarily  related  to  the  simplicity  of  the  mechanical  design. 

•  Very  little  quantitative  information  is  available  concerning  the 
fluid  mechanics  of  the  Coanda  eductor.  The  needed  information 
will  have  to  be  developed. 

•  Of  the  alternate  mathematical  approaches  available,  the  one  based 
on  an  existing  fluid  mechanic*  code  is  the  most  desirable. 

•  This  effort  should  be  continued  to  develop  the  information  base 
needed  to  design  Coanda  eductors  for  underground  excavation 
applications. 
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LITERATURE  REVIEW 


Literature  on  the  Coanda  effect  was  reviewed  with  the  purpose  of 
gathering  information  which  would  support  the  present  Coanda  eductor 
study.  A  bibliography  of  the  literature  reviewed  is  included  as 
Appendix  B.  The  review  revealed  that  very  little  information  is  avail¬ 
able  concerning  wall  jet  behavior  in  geometries  with  cylindrical  symme¬ 
try.  A  number  of  studies,  both  theoretical  and  experimental,  have  dealt 

with  two-dimensional  flow  over  flat  and  curved  surfaces.  In  addition,  a 

1 

relatively  large  body  of  information  has  been  developed  for  ejectors 
(or  jet  pumps)  which  entrain  fluid  by  flowing  primary  air  through  a  tube 
centered  in  a  laryer  mixing  tube.  For  the  Coanda  eductor,  the  primary 
air  is  introduced  through  an  annular  slit,  and  becomes  attached  to  the 
inlet  surface  prior  to  entry  into  a  mixing  zone.  Geometries  for  the 
two  eductors  are  similar  downstream  from  the  entry  to  mixing  section, 
hence  information  developed  for  this  region  in  ejectors  should  apply  to 
the  Coanda  eductor.  The  following  paragraphs  briefly  describe  the  infor¬ 
mation  obtained  from  the  literature  search. 

Studies  of  ejectors  have  included  both  theoretical  and  experimental 
approaches.  Theoretical  studies  consist  of  the  one-dimensional  approach, 
with  details  of  the  fluid  mechanics  typically  being  neglected.  Reversible 
thermodynamics  have  been  assumed  along  with  equations  of  state  for  the 
fluids.  Thus,  the  theoretical  analyses  establish  an  upper  limit  to  per¬ 
formance,  but  shed  little  light  on  how  to  design  Coanda  eductors  to 
achieve  the  maximum  performance.  Experimental  studies  with  air  ejectors 
have  shown  optimum  designs.  The  length-to-diameter  ratio  for  the  mixing 
section  is  found  to  be  near  7.  The  optimum  angle  of  divergence  in  the 
outlet  diffuser  is  found  to  be  in  the  7°  to  10°  range.  Because  of  the 
similarity  in  geometry,  these  optimum  parameters  would  be  expected  to 
apply  to  the  Coanda  eductor. 

Attachment  of  air  jets  to  surfaces  has  been  the  subject  of  a  number 
of  experimental  studies.  Applications  have  been  directed  primarily 
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toward  fluidic  amplifiers  and  to  aircraft  propulsion.  Jet  attachment  has 
been  studied  as  a  function  of  jet  width,  step  height,  and  angle  of  inclina¬ 
tion  of  attachment  surface  compared  to  the  slit.  From  the  results  one  can 
predict  the  distance  of  attachment  as  a  function  of  these  conditions. 

Several  studies  have  focused  on  the  effect  of  curvature  on  detachment  and 
fluid  entrainment.  Entrainment  increases  with  increasing  curvature  where 
curvature  is  defined  as  the  reciprocal  of  the  radius  of  curvature  (1/R). 
From  these  results,  one  would  expect  that  the  greatest  entrainment  ratios 
could  be  obtained  by  using  large  curvature  upstream  from  the  constant  area 
mixing  section  of  the  Coanda  eductor.  On  the  other  hand,  the  largest 
pressures  could  likely  be  developed  using  inlet  surfaces  which  employed 
less  curvature. 

To  date,  the  most  realistic  theoretical  approaches  have  relied  on 
similarity  theory  to  simplify  the  flow  equations,  thus  permitting 
solution.  This  similarity  approach  appears  to  be  quite  successful  in 
predicting  the  flow  field  for  flat  and  curved  two-dimensional  flows. 
Similarity  has  not  been  applied  to  the  axi symmetric  case  of  interest  in 
the  Coanda  eductor. 

Several  examples  of  reported  results  are  briefly  discussed  here. 

Figure  1  shows  the  results  of  Bourque  and  Newman^1  for  reattachment  of  a 
jet  to  an  inclined  plane.  The  minimum  length  of  plate  required  to  cause 
spontaneous  attachment  is  shown  as  the  upper  curve.  The  maximum  angle  at 
which  spontaneous  attachment  occurs  is  approximately  50°.  Once  the  jet 
becomes  attached,  however,  a  hysteresis  effect  becomes  evident  since 
attachment  can  be  maintained  up  to  about  6o°. 

Experimental  data  on  entrainment  for  jet  flow  along  the  surface  of  a 
cylinder  were  reported  by  Sridhar  and  Tn.^  Their  results  are  shown  in 
Figure  2,  in  which  the  dimensionless  entrainment  velocity  is  shown  as  a 
function  of  curvature  of  the  cylindrical  surface.  For  convex  surfaces, 
entrainment  per  unit  length  increases  as  the  radius  of  curvature  decreases. 
For  convex  surfaces,  just  the  opposite  occurs. 
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FIGURE  1  Results  of  Bourque  and  Newman^)  for  Attachment 
of  Plane  Jet  to  Inclined  Plane 
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FIGURE  3.  Cross  Section  of  Cylindrical  Eductors 
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MATHFMATICAL  model  development 


The  ultimate  use  of  the  mathematical  model  will  be  to  predict  mate¬ 
rial  transport  rates  for  eductors  employing  the  Coanda  effect.  The 
material  transport  capabilities  of  the  eductor  are  expected  to  be 
related  closely  to  the  air  flow  rate  and  the  pressure  drop  developed 
across  the  eductors.  Thus,  the  mathematical  model  should  permit  predic¬ 
tion  of  these  two  parameters  as  a  function  of  design  elements  including 
fluid  properties,  primary  jet  velocity,  and  geometric  design  of  a 
specific  eductor. 

To  develop  a  theoretical  model  within  the  time  allotted,  we  identi¬ 
fied  three  possible  approaches.  The  first  approach  was  to  use  model 
theory  which  relies  on  dimensional  analysis  to  permit  scaling  of  data 
from  model  to  prototype.  The  second  approach  was  to  use  similarity 
theory.  This  involves  simplifying  assumptions  to  convert  boundary  layer 
partial  differential  equations  to  total  differential  equations  which  are 
easily  solved.  Similarity  theory  has  been  used  successfully  in  modeling 
wall  jet  flows. (3)  The  third  potential  theoretical  approach  involved 
adapting  an  existing  fluid  flow  computer  code  capable  of  numerically 
solving  the  Navier-Stokes  equations  for  axial  symmetric  flow  fields. 

Our  first  efforts  were  to  briefly  survey  these  theoretical  approacnes 
to  verify  the  plausibility  of  each  for  application  to  the  present  prob¬ 
lem.  The  second  step  effort  was  concentrated  on  the  use  of  an  existing 
numerical  computer  code  since  it  promised  the  most  accurate  prediction 
of  the  flow  field  and  also  because  it  requires  the  least  amount  of  exper¬ 
imental  data  as  input.  At  present  we  have  experienced  some  success  in 
applying  the  computer  program;  thus,  we  are  continuing  to  concentrate  on 
this  approach.  If  the  numerical  approach  continues  to  prove  promising, 
no  additional  effort  will  be  devoted  to  the  other  two  alternates.  Each 
of  the  three  alternative  mathematical  approaches  is  described  further  in 
the  following  sections  of  this  report. 
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DIFFERENTIAL  EQUATIONS  OF  FLUID  FLOW 


The  mathematical  formulation  of  the  differential  equations  for  fluid 
flow  are  obtained  by  making  mass,  force,  and  energy  balances  on  a  small 
fluid  element  In  space.  These  basic  equations  form  the  basis  for  theore¬ 
tical  attacks  on  the  fluid  flow  problem.  Although  the  full  set  of  equa¬ 
tions  is  too  complicated  to  solve  generally,  they  were  listed  to  establish 
the  following: 

•  Simplifying  assumptions  used  in  various  theoretical  approaches 

•  Identification  of  important  physical  parameters 

•  Description  of  important  boundary  conditions 

•  Basis  for  dimensional  analysis. 

The  continuity  equation  results  from  a  statement  that  mass  Is  con¬ 
served.  For  cylindrical  coordinates,  tne  continuity  equation  Is 

It  +  r  &  (prvr)  *  F  -k  (pve)  *  £  h)  ■  0  o  > 

where 

p  =  fluid  density, 
t  *  time 

r  *  distance  measured  in  radial  direction, 
vf  *  fluid  velocity  in  radial  direction, 
ve  =  fluid  velocity  in  0  direction, 
vz  =  fluid  velocity  in  z  direction, 

0  *  angular  position  measured  around  z  axis, 
z  *  distance  measured  in  axial  direction. 

As  written.  Equation  (1)  applies  to  the  three-dimensional  flow  of  a 
compressible  fluid.  For  special  cases,  this  equation  may  be  appreciably 
simplified. 

The  momentum  equations  are  derived  by  balancing  the  forces  on  a 
fluid  element.  Statically  unbalanced  forces  are  equated  to  the  fluid  mass 
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sl'L?  aCCe,erati0"'  ThUS-  n  MJUnts  *»  application  of  Newton's 

the  n^rr:  ?  3  Sma11  f1Uid  e]ment‘  However*  deformation  of 
the  fluid  element  by  tangential  forces  must  also  be  accounted  for,  and 

nondefomaM  **  eqUat1°"S  °f  “>**>»  -pared  to  those  for 

Z  ™  In  LJt°r  rree;dt”e"Si0M'  f'°W-  three  fcaentum 

equations  result.  Written  In  cylindrical  coordinates.  In  terms  of  stress 

components,  the  momentum  equations  are  as  follows: 
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These  momentum  equations  apply  quite  generally  because  simplifying 
assumptions  have  not  yet  been  made. 
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The  energy  equation  Is  obtained  from  an  energy  balance  written  on  a 
small  fluid  element.  In  cylindrical  coordinates  It  may  be  written  In 
terms  of  stress  components  as 


where 

Cy  *  heat  capacity  at  constant  volume, 
T  =  temperature, 


%  = 


k2I 

3r  * 


%  = 


k  1  31 

r  30  ’ 


<*1 


»  • 


k  *  thermal  conductivity. 
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The  stress  components  which  are  included  in  Equations  (2)  through  (5) 
are  defined  as  follows. 


Equations  (1)  through  (5)  fully  describe  the  flow  of  fluids  in  three- 
space  dimension.  In  practice,  some  degree  of  simplification  is  required 
before  these  equations  can  be  solved. 
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DIMENSIONAL  ANALYSIS 


Any  mathematical  relationship  which  is  dimensionally  consistent  may 
be  written  as  a  relationship  between  dimensionless  groups.  The  various 
techniques  of  dimensional  analysis  identify  the  dimensionless  groups  for  a 
given  physical  problem.  The  utility  of  this  approach  is  that  it  can 
greatly  reduce  the  number  of  independent  variables  which  must  be  considered 
in  experiments.  Once  a  relationship  between  dimensionless  groups  is 
defined  by  means  of  model  experiments,  predictions  can  be  made  for  proto¬ 
types  provided  the  numerical  values  taken  by  the  various  dimensionless 
groups  fall  within  the  range  investigated  in  model  tests. 

One  way  of  deriving  the  Important  dimensionless  groups  for  a  problem 
is  based  on  simple  mathematical  operations  on  the  basic  differential  equa¬ 
tions.  This  method  is  well  described  by  Klinkenberg  and  Mooy.(4)  We  have 
used  this  approach  with  the  momentum  and  energy  equations.  Equations  (2) 
through  (5).  The  following  dimensionless  groups  shown  in  Table  1  appeared 
from  this  analysis. 

TABLE  1 .  Dimensionless  Groups  Obtained  from  Momentum 


and  Energy  Equations 

Name  of  Dimensionless  Group  Definition 
Reynolds  No.,  Re  pdv/y 

Prandtl  No.,  Pr  Cpy/k 

Eckert  No.,  Ec  vfycvt 

Euler  No.,  Eu  P/pv2 

Brinkman  No.,  Br  yv^/kT 

Specific  Heat  Ratio,  y  Cp/Cv 

Froude  No.,  Fr  9d/V^ 
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where 


p  =  fluid  density, 

V  *  velocity 
y  =  viscosity 

Cp  •  heat  capacity  at  constant  pressure 
k  *  thermal  conductivity 
Cy  =  heat  capacity  at  constant  volume 
T  =  temperature 
P  =  pressure 

g  =  acceleration  due  to  gravity 
d  *  diameter 

A  second  approach  to  deriving  important  dimensionless  groups  is  to 
list  all  important  parameters  for  a  given  physical  situation.  The  param¬ 
eters  can  then  be  grouped  as  dimensionless  ratios  using  standard  tech¬ 
niques.  A  hypothetical  Coanda  eductor  is  shown  schematically  in  Figure  5. 
In  addition  to  the  12  geometric  dimensions,  flow  properties  expected  ta  be 
important  include:  fluid  density,  fluid  viscosity,  specific  heat,  specific 
heat  ratio,  thermal  conductivity,  temperature,  velocity  of  primary  air, 
velocity  of  entrained  air,  and  pressure  drop  across  nozzle.  Using  the 
Buckingham  method^  for  forming  the  dimensionless  group,  we  obtain  the 
following  expression  for  the  entrained  air  flow. 

Entrained  Flow  ^  function  (Re,  Ma,  Br,  Eu,  y,  geometric  length  ratios)  (6) 
Primary  Flow 

where 

Re  *  Reynolds  No. 

Ma  =  Mach  No. 
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w.|  =  width  of  primary  air  inlet  slit 

w2  =  height  of  step  separating  jet  and  Coanda  wall 

w^  =  length  of  inlet  cylinder  upstream  from  primary  air  slit 

<*1  =  angle  of  divergence  of  outlet  diffuser 

a2  =  angle  of  inclination  of  primary  jet 

a^  =  angle  of  convergence  of  inlet 

=  length  of  mixing  section  (constant  area) 

£2  =  length  of  outlet  diffuser 
£3  =  length  of  Coanda  entrainment  section 
£4  =  length  of  inlet  transition 
Rl  =  radius  of  throat  of  nozzle 

R2  =  radius  of  curvature  of  Coanda  entrainment  section 

FIGURE  5.  Schematic  View  of  Hypothetical  Coanda  Eductor 
Showing  Geometric  Parameter 
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Br  ■  Brinkman  No. 

Eu  a  Euler  No. 

Y  *  specific  heat  ratio. 


The  Mach  No.  is  closely  related  to  the  Eckert  No.  which  appeared 
earlier.  From  Equatloo  (6).  one  could.  In  theory,  perform  experiments ;  that 
would  provide  the  functional  relationship  between  the  groups,  n  pr 
;:7ny  experiments  wou'd  be  required  to  evaluate  the  relationship  among 
the  five  flow  related  groups  and  the  12  geometric  ratios.  The  next  step 
in  using  the  dimensional  analysis  approach  would  be  to  eliminate  the  g  P 
of  minor  Importance,  and  key  on  the  several  controlling  groups.  This  step 
was  not  accomplished  in  the  present  work. 

SIMILARITY  THEORY 

similarity  theory  Is  a  method  for  simplifying  the  equations  of  motion  to  a 

y  k  eniwpH  The  term  "similarity"  derives  from  the 

point  where  they  can  be  solved.  The  term  sum*  y 
assumed  shape  of  the  velocity  profile,  which  is  >«umed  to  ave 
scribed  similarity  at  each  downstream  position,  x.  A  c  as  leal  e P  _ 
s  mUarity  theory  is  that  used  by  Glauert'3'  In  his  analysis  of  t  e  wall 

.  t  ic  the  simDlified,  two-dimensional  boundary 
jet.  Glauert's  beginning  point  is  the  simpiiT 

layer  momentum  equation  written  as 


su  t  v  BL  . 

3X  3y  9y 


(7) 


where 


u  =  velocity  in  x  direction, 
v  =  velocity  in  y  direction, 
v  =  kinematic  viscosity. 
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A  stream  function,  *,  Is  defined  satisfy  the  continuity  equation: 


24 

9x 


XV 


and 


9  y 


XU 


(8) 


(9) 


where 


ip  =  stream  function. 


The  continuity  equation  may  be  written  as 


£x  (xu)  +  ^  (xv)  =  0  (10) 

By  defining  a  new  variable,  n  =  yx11,  the  two  partial  differential  equa¬ 
tions,  Equation  (7)  and  Equation  (10)  are  transformed  to  a  total  differen¬ 
tial  equation: 

f‘"  +  ff"  +of'2  =  0  (1]) 

In  Equation  (11)  the  argument  of  the  function,  f.  Is  n;  a  Is  a  numerical 
parameter  related  to  n. 

(3 ) 

Glauert's  solution  to  Equation  (11)  provides  velocity  profiles 
which  agree  well  with  experimental  measurements,  showing  that  the  mathe¬ 
matical  assumptions  made  were  justified.  Presumably,  similar  approaches 
could  be  made  for  ax i symmetric  flows  encountered  in  the  Coanda  eductor. 

No  attempt  was  made  to  pursue  the  similarity  theory  approach,  but  it 

appears  to  offer  potential  and  serves  as  a  backup  in  the  event  the  numeri¬ 
cal  approach  should  fail. 
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NUMERICAL  MODEL  BASED  ON  VORT1CITV  TRANSPORT  THEORY 


The  following  discusses  the  mathematical  basis  and  the  i'nite 
difference  method  of  solution  used  in  the  vorticity  transport  method. 

Mathematical  Basis 

The  differential  equations  governing  Coanda  eductor  fluid  flow  are 
derived  from  the  two  conservation  laws: 

•  conservation  of  mass  (continuity) 

•  conservation  of  momentum  (Newton's  second  law). 

These  equations  written  for  incompressible,  turbulent  flow  in 
cartesian  tensor  form  are: 

Continuity: 


Momentum: 


(12) 


1  JE_ 
p  ax. 


(13) 


In  the  above  equations,  time  and  spatial  changes  of  the  fluid  density, 
p,  have  been  ignored  (incompressibility  condition). 


Notation  used  in  Equations  (12)  amd  (13)  is  as  follows: 

i.  n 

u.j  =  velocity  along  itn  coordinate 

xi  =  itn  space  coordinate 

P  =  pressure 
t  =  time 

e..  =  Eddy  diffusivity  tensor  for  momentum. 
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The  operator  D/Dt  is  the  substantial  derivative  defined  as 


_D_  _  _3_  a 
Dt  “  at  ax. 


The  tensor  index  i  (or  j)  for  three-space  takes  values  i  =  1,  2  and  3. 
Einsteinian  notation  is  used  where  a  repeated  index  implies  summation  over 
the  possible  values  i  (or  j)  =  1,  2  and  3. 

In  the  momentum  Equation  (13),  we  have  neglected  gravitational  effects 
and  have  assumed  that  the  effects  of  Reynolds  stresses  (turbulent  stresses) 
may  be  approximated  through  the  use  of  an  eddy  diffusivity,  e. .. 

•  J 

For  the  Coanda  eductor  study,  it  is  possible  to  express  the  equations 
of  continuity  and  momentum  in  terms  of  axisymmetric  coordinates,  thus 
using  only  the  two-space  coordinates  z  (axial)  and  r  (radial).  This 
assumption  ignores  mean  azimuthal  variation  of  the  dependent  variables  such 
as  mean  azimuthal  velocity  (swirl).  Also,  the  Coanda  study  is  directed  to 
steady- flow  operation  of  the  eductor,  hence  time  dependence  may  also  be 
eliminated  in  the  equations  of  motion. 

Thus,  the  equations  which  are  used  for  analysis  of  the  Coanda  eductor 
are  the  following  steady-flow,  axisymmetric  equations  of  continuity  and 
motion: 

Continuity: 


1  lUL  +  3V. 
r  3r  az 


=  0  . 


(14) 


Momentum: 

r-direction, 
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z-direction. 


In  the  above  equations,  u  and  v  are  velocity  in  the  r  and  z  direc¬ 
tions,  respectively,  and  the  eddy  diffusivity  tensor  is  assumed  to  have 
components  er  and  ez  only.  With  the  present  state-of-the-art,  it  is 
virtually  impossible  to  treat  the  diffusivities  with  more  detail. 

Use  of  only  the  er  and  ez  components  of  momentum  diffusion  is  to 
assume  that  the  diffusion  coefficient  for  v-velocity  and  u-velocity  are 
the  same  in  a  given  direction. 

Vorticity  Transport 

In  the  application  of  numerical  analysis  we  consider  solving  Equa¬ 
tions  (10.  (15)  and  (16)  by  finite-difference  techniques  without  further 
simplification.  However,  we  must  devise  some  method  to  solve  for  pres¬ 
sure,  P.  This  may  be  done  by  taking  the  divergence  of  Equations  (14)  and 
(15)  and  adding  the  results  to  obtain 


where  the  operator 


32  .  1  9  .  32 

2  2 
3r  r  3r  3z 


In  Equation  (17),  terms  involving  derivatives  of  er  and  ez  have  been 
ignored  for  the  sake  of  this  discussion.  To  obtain  a  solution  to  the 
Coanda  dynamics.  Equations  (15),  (16)  and  (17)  would  be  solveo  simulta¬ 
neously.  However,  experience  has  shown  that  Equation  (17)  is  difficult  to 
solve  numerically  because  of  required  boundary  conditions. 
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To  circumvent  simultaneous  solutions  to  the  pressure  equation  and 
associated  difficulty  with  the  boundary  condition,  one  may  choose  to  solve 
an  equivalent  set  of  equations  which  eliminates  pressure  as  an  explicit 
part  of  the  analysis.  This  approach  is  called  the  vorticity  transport 
method . 


We  define  vorticity,  id,  as 


co 


3u  3v 
3z  "  3r  * 


A  stream  function,  if >,  is  defined  by 


“  ‘r  It  arK* 

r  3r 


1M 

r  3z 


(18) 


(19) 

(20) 


which  identically  satisfies  the  continuity  Equation  (14).  An  elliptic 
partial  differential  equation  is  obtained  for  \p  by  substituting  Equa¬ 
tions  (19)  and  (20)  into  (18)  which  yields 

4-  mtik.  .ru.  (21) 

3r  r  3r  3z 

Velocities  u  and  v  may  be  obtained  by  solving  Equation  (21)  for  i p  and  then 
solving  the  auxiliary  Equations  (19)  and  (20). 

However,  to  obtain  solution  to  Equation  (21),  the  flow  field  vorticity, 
w(r,z),  must  be  known.  A  vorticity  transport  equation  may  be  derived  by 
cross-differentiating  Equations  (15)  and  (16)  and  subtracting  the  latter 
result  from  the  former  to  obtain 
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(22) 


If  the  turbulent  structure  of  the  flow  f’eld  is  homogeneous  and  iso¬ 
tropic,  derivatives  of  er  and  ez  vanish  and  the  vorticity  transport  equa¬ 
tion  becomes 


3uu  3vu 
3r  3z 


3  /l  3a)r\  3  oj 

*<>  —  \~r~  +E,=  r7 


3r  \r  3r 


3z 


(23) 


Note  that  pressure,  P,  does  not  explicitly  appear  in  Equations  (22)  and 
(23).  Once  the  vorticity  solution  is  obtained  (and  hence,  the  velocity 
solution),  one  may  back-calculate  pressure  from  Equation  (17).  Alterna¬ 
tively,  pressure  may  also  be  calculated  using  either  Equation  (15)  or  (16). 

Description  of  the  SYMJET  Code  (Coanda  Version! 

The  Coanda  version  of  the  SYMJET  code  solves  Equations  (21)  and  (22) 
along  with  the  auxiliary  relationships  Equations  (19)  and  (20),  by  finite- 
difference  techniques.  However,  these  equations  are  first  scaled  with  the 
dimensionless  variables  and  parameters  defined  as  follows: 

R  =  r/rQ 

Z  =  z/rQ  ,  I  =  z/D 
U  =  u/vQ 
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V  =  v/v0 

*  ‘  */(r0  vo) 

1  =  b)  r0/v0 

P  V 

RE„  =  ,  (radial,  turbulent  Reynolds  No.) 

r  er 

X*  V 

RE,  =  ,  (vertical,  turbulent  Reynolds  No.) 

z  e2 

In  the  above,  rQ  is  the  characteristic  radius  of  the  Coanda  eductor  and  v( 
is  a  characteristic  reference  velocity  which  may  be  taken  as  the  Coanda 
slot  velocity. 

In  their  scaled  form  the  governing  equations  become: 

Stream  Function: 


1  av 


2 

31 


R  SR  3? 


Rn 


Vorticity  Transport: 


aun  avn  i  / 32c  i  an 

—  +  —  =  —  —j  + - 

aR  az  RE,,  R  3R 

r  \ 
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2 

a^n 


rez  az2 


along  with  the  auxiliary  relationships  for  velocity: 


U  = 


1  ai 

r  az 


v 


1 II 
r  aR 


(24) 


(25) 


(26) 

(27) 
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Terms  involving  derivatives  of  REr  and  REZ  are  not  included  in  Equa¬ 
tion  (25)  but  are  accounted  for  in  the  SYMJET  code. 

Finite  Difference  Grid  System 

The  finite  difference  grid  layout  consists  of  two  grid  systems.  One 
grid  is  used  to  calculate  the  stream  function,  ¥,  which  provides  informa¬ 
tion  to  compute  velocity  components,  U  and  V.  This  system  coincides  with 
the  physical  bouidaries  and  is  illustrated  by  the  wider  lines  on  Figure  6. 
The  stream  function  is  calculated  at  the  system  interior  intersection 
points  designated  by  the  solid  round  symbols.  The  U  components  of  the 
velocity  field  are  computed  at  vertical  midpoints  which  are  designated  by 
open  circle  symbols;  whereas,  the  V  components  are  computed  at  horizontal 
midpoints  (K  coordinate)  and  designated  by  open  box  symbols.  In  this 
manner,  the  stream  function  grid  layout  defines  a  system  of  cells  with  the 
stream  function,  ¥,  computed  at  each  corner  point  (or  set  by  boundary 
conditions,  as  the  case  may  be)  and  velocities  defined  at  the  center  of 
the  cell  face. 

The  second  griu  system  is  used  to  calculate  vorticity,  n,  and  is 
illustrated  in  Figure  6  by  the  narrow  lines.  This  layout  completely  over¬ 
laps  the  ¥  grid  (and  physical  system)  with  interior  intersection  points 
centered  in  the  cells  which  are  defined  by  the  ¥  grid  system.  These 
interior  grid  points  are  indicated  by  crosses. 

This  staggered  grid  system  is  used  for  computational  convenience  in 
treating  boundary  conditions  and  to  permit  direct  evaluation  of  convective 
transport  terms  at  cell  faces. 

The  ¥  grid  system  is  sized  by  NJ  and  NK  grid  points  in  the  R  direction 
and  vertical  direction,  respectively  The  n  system  has  size  NJ  +  1  and 
NK  +  1  in  the  respective  directions.  Points  on  the  ¥  grid  are  indicated  by 
(ji  k),  whereas  points  on  the  n  grid  are  specified  by  (p,  q).  Vertical 
spacing  for  the  system  is  defined  by  AZk  which  may  be  variable.  Grid 
spacing  along  the  R  coordinate  is  designated  by  AR. 
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Difference  Equations 

Standard  difference  representation  is  used  wherever  possible  in  this 
discussion.  Central  differences  are  used  for  both  first  and  second  partial 
derivatives  except  for  convective  terms  where  special  upstream  methods  are 
used.  These  methods  are  similar  to  those  used  by  Torrance  and  Rockett^ 
and  Runchal  and  Wolfshteln,^  and  details  of  the  application  here  are 
explained  by  Trent  ^  Techniques  for  uneven  spacing  are  used  for  the 
vertical  differences. 

Stream  Function  and  Velocity 

Consider  the  stream  function  grid  system  illustrated  ir.  Figure  6. 

The  finite  difference  representation  of  Equation  (24)  based  on  central 
differences  for  both  first  and  second  partial  derivatives  is  as  follows: 


2 


+ 


+ 


TJ.k-l 


(28) 


In  the  above  difference  equation,  the  quantity  ^  is  the  average  value 

vorticity  at  point  (j,k),  hence  the  overbar.  This  average  value  is  used 

since  does  not  lie  on  the  f  computational  grid. 

P»9 

Velocity  is  calculated  at  cell  faces  by 
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(29) 
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and 


Transport  Equations 

Referring  to  the  p.q  grid  system  Illustrated  In  Figure  6.  the  differ¬ 
ence  representation  of  the  steady  flow  vortlclty  transport  Equation  (25)  is 
written  as  (after  collecting  terms). 
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The  SYMJET  coda  solves  the  two  difference  Equations  (28)  and  (31)  by 
Gauss-Siedel  iteration.  Iteration  Is  continued  until  the  computed  velocity 
meets  a  particular  convergence  criterion.  Although  this  criterion  is  left 
to  the  discretion  of  the  user,  computation  in  the  present  work  was  carried 
out  until  the  maximum  relative  change  in  velocity  for  the  entire  flow  field 

_3 

was  less  than  10  for  one  Iteration. 

Application  of  SYMJET  Code  to  Coanda  Eductor 

The  SYMJET  computer  code  is  presently  being  applied  to  predict 
detailed  flow  fields  in  the  Coanda  eductor.  Initial  efforts  in  applying 
the  code  were  to  develop  a  routine  to  permit  arbitrary  specification  of 
boundary  conditions  from  input.  This  routine  has  proven  highly  successful. 
At  present,  all  boundary  conditions  can  be  easily  specified.  A  search 
routine  built  into  the  program  maps  the  specified  points  onto  the  overall 
grid  system. 

At  present,  some  success  has  been  achieved  in  predicting  flow  fields 
for  realistic  cases.  Steady  state,  well -converged  solutions  have  been 
obtained  for  several  flow  situations.  The  grid  system  used  for  these  cases 
is  showc  in  Figure  7.  The  solid  boundaries  were  selected  to  closely  model 
an  existing,  workable,  Coanda  eductor. 

At  the  center  of  the  nozzle,  the  boundary  condition  requires  that  the 
velocity  gradient  in  the  radial  direction  be  zero.  Along  the  outer  walls 
of  the  eductor,  the  no-slip  criteria  (velocity  zero  at  olid  wall)  was 
used.  At  the  outlet,  velocity  profiles  were  specified;  at  the  inlet  (top 
of  Figure  7),  the  free-flow  boundary  condition  was  used.  This  condition 
requires  only  that  the  streamlines  remain  straight  as  they  enter.  At  the 
primary  jet  opening,  the  velocity  was  specified. 

The  unequal  grid  spacing  evident  in  Figure  7  was  required  to  ade¬ 
quately  handle  flow  issuing  from  the  primary  jet.  Initial  attempts  to  use 
an  equal  spacing,  with  only  a  single-mesh  node  being  within  the  primary 
jet,  failed  to  converge. 

Only  the  inlet  section  of  the  Coanda  eductor  was  included  in  the  cal- 
culational  grid  system.  Although  it  is  possible  to  include  downstream 
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FIGURE  7. 


Calculatlonal  Grid  System  Showing  Boundary 
Conditions  and  Grid  Spacing 
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regions,  this  would  substantially  increase  the  number  of  calculatijnal 

nodes  and  the  computer  time.  This  was  considered  unnecessary  for  initial 
runs. 


Results  for  a  typical  case  in  which  the  total  air  flow  was  five  times 
the  flow  of  the  primary  air  are  presented  in  Figure  8.  The  shape  of  the 
streamlines  demonstrates  that  the  primary  jet  becomes  attached  to  the 
outer  wall  of  the  eductor,  and  entrains  additional  air. 

Figure  9  shows  results  obtained  for  a  lower  entrainment  ratio.  The 
total  air  flow  was  reduced  by  lowering  the  specified  outlet  velocity.  The 
total  flow  leaving  through  the  outlet  (bottom  of  Figure  9)  was  slightly 
lower  than  that  entering  through  the  primary  air  slit.  Flow  separation 
occurred  near  the  center  of  the  nozzle,  with  a  part  of  the  primary  air 
exiting  through  the  normal  Inlet.  It  Is  worth  noting  that  this  type  of 
behavior  could  not  be  predicted  using  similarity  theory  because  the  /elocity 
profiles  are  not  similar  at  the  various  downstream  positions. 

In  addition  to  the  streamline  output  shown  in  Figures  8  and  9,  the 
SYMJET  program  prints  velocity  components,  stream  functions,  and  vorticity 
at  each  flow  node. 

For  the  cases  shown  In  Figures  8  and  9,  the  eddy  viscosity  was  taken 
as  constant  throughout  the  whole  flow  field.  This  is  recognized  as  a 
gross  simplification  which  will  be  Improved  upon  as  soon  as  comparative 
experimental  data  become  available. 

Pressure  is  eliminated  as  a  variable  In  the  vorticity  transport  method 
of  flow  calculation,  hence  It  was  not  calculated  In  predictions  made  to 
date.  Since  pressure  drop  Is  an  Important  consideration  in  the  operation 
of  an  eductor,  the  next  major  effort  in  this  work  will  be  devoted  to 
numerical  calculation  pressure. 
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Streamlines  Predicted  for  an  Entrainment  Ratio  of  0.8 
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EXPERIMENTAL  PROGRAM 


The  experimental  apparatus  is  shown  in  Figure  10.  It  will  be  used  to 
gather  data  on  pressure,  temperature,  air  flow  and  velocity  profiles  at 
critical  sections  of  the  Coanda  eductor.  Performance  predictions  of  the 
analytical  program  which  will  be  matched  by  the  experimental  program  are: 
tr.trained  flow  rate  to  primary  jet  flow  rate,  maximum  total  developed  pres¬ 
sure,  static  pressure  difference  between  eductor  inlet  and  outlet,  velocity 
profiles  in  the  entrainment  sections  of  the  eductor,  and  efficiency  of 
energy  Interchange  of  the  primary  air  jet  to  combined  primary-entrained  air 
stream. 

The  initial  experiments  required  in  support  of  the  analytical  model 
will  provide  data  that  will  allow  the  determination  of  turbulent  viscosi¬ 
ties  in  the  inlet  region  of  a  Coanda  eductor.  This  Information  is  essen¬ 
tial  to  obtain  realistic  predicts  nr,  because  turbulence  is  expected  to  play 
an  Important  role  in  attachment  and  detachment  of  the  primary  jet  and  to 
entrainment  of  fluid. 
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hematic  View  of  Experimental  Equipment 
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